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1.0 


INTRODUCTION AND SUMMARY 


Two techniques for simultaneously estimating altitude, ocean 
wave helgnt, and signal-to-noise ratio from the 6E0S-C satellite 
altimeter data are described In this report. One technique Is based 
on maximum likelihood estimation, MLE, and the other on minimum 
mean square error estimation, MMSE. Performance is determined by 
comparing the variance and bias of each technique with the variance 
and bias of the smoothed output from the GEOS altimeter tracker. 

Ocean wave height tracking performance for the MLE and MMSE algorithms 
Is measured by comparing the variance and bias of the wave height 
estimates with that of the expression for t^e return waveform 
obtained by a fit to the average output of the 16 waveform sampling 
gates. 

Both the MLE and MMSE algorithms achieve Improved altitude 
tracking performance over the GEOS trackers with real data In 
simulated "real time". In fact, as shown in the results (see 
Section 2.0), both algorithms achieve close to theoretical performance 
in all parameter estimates (altitude, wave height, signal and noise). 

It Is significant to note that both algorithms out-perform the GEOS 
tracker despite the fact that the processing was performed on a 
poor data tape with low sea state and only the 16 leading edge 
samples could be used to provide the parameter estimates. Thus, 
it Is felt that with good data from an altimeter designed to accommodate 
the MLE and MMSE algorithms (better range resolution and more leading edge 
samples as well as usable noise and plateau gate samples), the 
algorithm could achieve the theoretically predicted factor of four 
reduction in the standard deviation of the altitude estimated as 
compared to the standard split-gate tracker at high sea states. 
Furthermore, with such a properly aesigned altimeter, the standard 
deviation of the other parameter estimates (wave height, signal-to-noise 
ratio, and pointing error) could be expected to be on the order of 
1% of the parameters me*r> value. 


A major problem was encountered in using the MLE and MMSE 


algorithms for ground processing the GEOS altimeter data. Since 
In GEOS, the locatlo In range of the 16 sample gates Is varied 
according to the detected altitude error on a per pulse basis . It Is 
necessary to process via the MLE or MMSE algorithms on a per pulse 
basis. Thus It would be Impractical, from a computational standpoint, 
to utilize either the MLE or MMSE algorithms as an operational GEOS 
ground data processing system. For this reason, the GEOS data runs 
that have been processed have been very short and are Intended to 
verify the algorithms performance. However, In future altimeters 
the range samples can be preaveraged and the algorithms described 
In this report could be readily Implemented as either an on-board 
data processor or as a general data processing system. 
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2.0 


RESULTS OF GEOS-C DATA PROCESSING 


The performance of the MLE and MMSE algorithms is summarized 
In Table 1 for two different portions of orbit 186. The GEOS 
split-gate tracker performance Is also shown to provide a direct 
comparison of the standard deviations in the altitude estimates. As 
shown in the table, the improved tracking algorithms reduce the standard 
deviation in the altitude estimate by about 25 to 30 centimeters, or 
30 percent, compared to GEOS. While this In Itself Is certainly 
not considered a significant improvement, it should be noted that 
the measured performance is very close to the theoretically predicted 
performance. Thus, while the GEOS altimeter design Itself limits 
the Improved performance, achieving near theoretical performance 
with a properly designed altimeter would result In a factor of four 
reduction in standard deviation at high sea states. Of course, 
considerably more data runs (especially at different sea states) 
would be required to ultimately verify this conclusion. 

The following comments are appropriate for a more comprehensive 
understanding of the data presented in Table 1. 

1. The number of pulses averaged in the tracking loops is 
7.56 for altitude (MLE and MMSE matched to GEOS) and 40.8 for the 
other three loops in all cases. 

2. The results of two runs are shown; one containing files 
69-72 and the other files 10-14. Each runs contains 12 seconds 
of data (1200 pulses). 

3. For each run, a least squares straight line fit is made 
to the GEOS instantaneous altitude data. The measured altitude 
standard deviations on a per pulse basis for GEOS, MLE and MMSE 
are then computed from the residues about this line. The line is 
also used to determine the biases for the MLE and MMSE relative to 
GEOS. As an example, the straight line for files 69-72 is 

y 3 4.46 + 13.5155 AT meters where AT is the interpulse period and 
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RESULTS OF GEOS-C DATA PROCESSING FOR ORBIT 186 
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NOTE: * Values in italics correspond to a static MLE fit to the data and are used to calculate 
theoretical performance. 






















the line is relative to an altitude of 853,300 meters. 


4. The MLE and MMSE theoretical performance is based on a 
numerical evaluation of the covariance matrix of the parameters. 

To perform the numerical evaluation. It Is assumed that the static 
fit parameter values are the correct values, and that the 16 range 
samples are independent. The GEOS theoretical performance Is based 
on the previously derived closed form expression for split-gate 
trackers and the static fit parameter values. For GEOS, It Is also 
assumed that there are two Independent samples In both the early 
and late gates, which is consistent with the 16 Independent sample 
assumption for the MLE and MMSE. 

The need for additional data processing at other sea states 
Is readily demonstrated in Figure 1, which shows the measured 
standard deviation in the altitude estimate versus measured wave 
height for the MLE and MMSE and versus the static fit wave height 
for GEOS. 

To achieve some degree of confidence in the measured mean parameter 
estimates (as well as for use in evaluating the theoretically 
predicted standard deviations), it was decided to perform a static 
MLE fit to one of the data runs (files 69-72). This was done as 
follows: The GEOS instantaneous (per pulse) altitude estimates 

were quantized into 23 range bins, each of 1/4 gate width (23.44 cm), 
centered about the mean straight line fit to GEOS altitude. The 
resulting distribution of C.te #10 is shown In Figure 2. Now, the 
resulting returns (each of the 16 range samples) for each bin were 
averaged and then a MLL performed on this averaged return for each bin. 
Figure 3 shows the MLE fit for bin zero. Finally, a weighted 
(proportional to number of pulses in a bin) average of the 23 MLE 
estimates is computed to obtain what is referred to as a "static 
MLE fit" to the data. These are then the parameter values shown in 
Table 1. Considering the excellent agreement between the static fit 
values and the measured values and the effort required to perform a static 
fit to the data, only measured parameter values should be used in any 
future data runs. 
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3.0 


DEVELOPMENT OF THE MLE AND MMSE ALGORITHMS 


Consider the altimeter receiver which can be modeled as (nearly) 
matched filtering, square law detection and sampling as shown in 
Figure 4. The description of maximum likelihood estimation* 
and Its properties are simplified with this model, since the problem Is 
now constrained to that of finding the optimum estimates for the 
sampled video outputs. All that Is required Is a parametric expression 
for the mean video power return, V. In general, V, is a function 
of a large number of parameters. These Include not only the usual 
radar parameters required to establish receiver slgnal-to-nolse ratio, 
but also antenna beamwldth, off-nadir pointing error, range resolution, 
and the height and slope distributions of the ocean surface. The 
trick Is to reduce the number of parameters required to specify 7 
to a minimum, since a joint maximum likelihood estimate of all of 
these parameters will be necessary. It can be shown that an 
approximate closed form solution for V can be developed that Includes 
all of the above parameters. That Is, the average video power return 
versus range y from the altimeter can be expressed In the fonr. 


where 


and 


V(yly 0 .a,o e ,u) (1) 

I 

y Q - Is the altitude (range) of the a tlmeter above 
the mean sea surface; 

a - Is the slgnal-to-nolse r»tio; 

o p - Is the effective RMS waveheight; 

p - Is the exponential decay factor. 


Later It will be shown that the MMSE estimation algorithm Is the 
same as the MLE except that the return difference Is weighted by 
unity Instead of Its variance. Thus the following discussion Is 
appropriate to both the MLE and MMSE algorithms. 
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Explicit expressions for a, o g and p 1r terms of physical parameters are 
given In Appendix A. 


3.1 Estimation Equations 

The procedure to be followed In developing the maximum likelihood 
estimation equations will be to find a suitable approximation tc 
the likelihood function (l.e., the joint probability density) o f 
the video outputs In terms cf the four parameters y Q , a, and p. 

Then for any set of observed video samples, toe joint maximum likelihood 
estimates are the values (y Q , a, o e , p) which maximize the logarithm 
of the likell ood function. The maximum, of course, occurs at 
a point 1 r. the parameter space for which the four partial derivatives 
with respect to the parameters are zero. 

The form of the likelihood function for the sampled video 
outputs of Figure 1 Is simple, if It Is assumed that the samples 
are separated by at least one range resolution cell. In this case, 
the correlation between samples Is negligible, and since the underlying 
process Is Gaussian, It Is reasonable to assume them to be Independent. 
Thus, the likelihood function has the form: 

-1 

A(V 11 ,V 12 ,...V nk /y 0 ,a,o e ,p) * 1 n |( V 1 exp (-V 1k /V <| ) (2) 


where: 

V . is the sampled video output from the 1— 

IK 4L. 

range cell on the k-- pulse. 


and 


V 1 is the expected value of the video from the 1— range cell. 

That is, the square law detected video outputs form an independent 
exponential process. The dependence of the likelihood function on the 
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four parameters y Q , a, o e , u Is contained entirely In the vari- 
ation with range of the mean value Vj . Thus taking the logarithm 
of (2), and differentiating, the maximum likelihood estimates 
must satisfy the four equations 


0 . AlfiSA --is 

k 1 


- -2 




< v i - v i k > v t 


where a stands for any of the four parameters y Q , a, o g , p. 


( 3 ) 


3.2 Im plementation 

The problem of determining the joint maximum likelihood 
estimates of the four parameters describing the altimeter video 
power return requires solving Equations (3). One technique for 
achieving this Is to use the negative of the partial derivatives 
in Equations (3) as inputs to Integrating - feedback filters to 
derive the estimates as shown In Figure 5. 

That is, since the estimates are derived from Integrating 
filters, one has: 


6 ls. . a'QQA 

3T ,y 0 


da 

Ht 


ja 


d<je _ . 3 log A 

dt 35 

e 


(4) 


where k is a gain factor associated with the integrators. Now, the time 
rate of change of the logarithm of the likelihood function is: 
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Figure 5. General MLP Implementation 
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( 5 ) 


< 0 . 

Thus, as a function of time, the logarithm of the likelihood 
ratio must decrease and converge to a local minimum of the likelihood 
function. If tne initial parameter estimates are close to the 
correct ones, then the loops In Figure 5 will converge to the joint 
maximum likelihood estimate of y Q , a, o p and u. 


The practicality of implementing the MLP depends on the 
complexity of the partial derivative calculator In Figure 5. 

The steps required to compute the four error signals are diagrammed 
In Figure 6. As can be seen from that figure, if the mean video 
versus range plus the four partial derivatives are available, 
then the computations required to derive the error signals are 
relatively minor. 


Figures 5 and 6 give a fairly pleasing intuitive picture of 
the operation of the MLP. In Figure 6, it is seen that the error 
signals for the tracking loops in Figure 5 are derived as follows: 
First, an error signal versus range, V - V, is generated which 
represents the difference between the measured return as a function 
of range, V, and the estimated mean return, V. Secondly, this error 
signal is weighted inversely proportional to its variance*. That 
is: 


E ( V - V) 2 = V 2 (6) 

Thus, after normal ization, the signal may be loosely described as having 
uniform information content. Finally, the normalized error-versus- 
range signal is "gated" by the partial derivative of V, with respect to 


*For MMSE estimation, the error signal weighting is unity instead 
of its variance. i/i 


f^LLZAU 



Figure 












the parameter of Interest, and summed over range. The effect of 
multiplying (gating) by the partial derivative Is to emphasize 
those range bins which are most affected by variations In the 
parameter of Interest. 


3.3 Decoupling Matrix 

A major shortcoming of the feedback loop Implementation described 
above Is that the loop error signals are very highly cross-coupled. 

As a result, an error in one of the parameter estimates can (and 
usually does) produce errors In the other estimates. Furthermore, 
experience with the Implementation has shown the convergence time 
of the loops to be excessively long as a result of the error 
cross-coupling. These difficulties may be essentially eliminated 
by decoupling the loop error signals. Fortunately, a decoupling 
matrix can be Implemented with a fairly simple technique. 

The underlying assumption In the decoupling matrix approach 
is that the loop error signals and the true errors are linearly 
related. However, since the loops are very highly cross-coupled, 
the linear relation must be described by a full matrix transform: 

£ * KAa (7) 

where c - Is a 4 x 1 vector of error signals which are generated 
by the true errors Aa(Au is a 4 x 1 vector). 

and K - Is a 4 x 4 matrix describing the linear relation. 

The basic problem in determining the matrix K Is that it depends on 

the (4 x 1 vector) of true parameter values, a s La, y Q , o e » u] T which are, 

of course, unknown. An obvious approach to solving this problem would 

be to evaluate the matrix K at the estimated parameter values, a and 

use the estimated matrix to decouple the loops. The usefulness of 

this approach will, of course, depend on the region over which a 

linear relation, Equation (7), holds. If the region is very small, 
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then a would have to be very close to the true value a. Here again, 
experience with the Implementation has shown that over an acceptably 
large region about the true values, a linear relation exists 
between an Input error In one parameter and the corresponding errors 
produced In the other three parameters. 

Next, consider the problem of estimating the matrix K. Suppose 
there Is an error In a single parameter, say Aw. If the error 
signals are linearly related to the true error, then: 



Where the notation e^ lJ Is used to represent the error signal produced 

a 

in the "a" loop by an error of Aw In the w parameter. From the 
above equation, It Is seen that the fourth column of the K matrix 
Is given by 
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Hence the K matrix can be computed according to the relation 


K - 

0/A«) 

r*!*" 

• • • (1/ A|i) 





m 

c 

« 






Ax 

« 0 


• 




e 

Aa 


A'i 




w » J 





- 






The procedure for decoupling the loops Is as follows. 

1) After the Initial large transients have died down (—4 iterations) 
compute a simulated return based on the current parameter estimates. 

2) Using the simulated return as the correct return, perturb the 
parameter estimates one at a time and compute the corresponding 
error signals. 

3) From these four sets of error signals, build up the K matrix 
according to Equation (10). 

4) Invert K, and there after decouple the error signals according 

to the formula: Aa * K c(o) 


In practice, the true parameter values will be changing slowly, so 
that this procedure must be repeated at periodic Intervals. Figure 7 
shows a functional block di ;ram of the resulting feedback loop 
implementation of the MLP when the error decoupling matrix is included. 


Theoretical Accuracies 


Finally, it is useful to derive the general expressions for the 
theoretical accuracies (asymptotic variance) of the joint maximum 
likelihood estimates. These accuracies are given in terms of 
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Figure 7. Functional Block Diagram of the MLP 









the standard deviation of the estimates due to random fluctuations 
only . System errors (such as clock and timing errors) as well as 
errors associated with spatial averaging are not Included In the 
evaluation. 

From the theory of maximum likelihood estimation (see for example 
Wilks [1], Chapter 12), It Is known that for regular distributions, 
the MLE is, asymptotically, unbiased and has minimum variance. 

The covariance matrix of the estimation errors Is given by 

C « B _1 (11) 


where C Is the covariance matrix defined by: 


C, m " Covar (a. « m ) 
4m *• m 


with <j,a m standing for any of the parameter estimates and B Is the 
Information matrix defined by: 



d log A 


d a. 


d log A 

da 

m 


( 12 ) 


where. In (12) A stands for the likelihood atio evaluated at the 
true parameter values. 


Substituting (3) into (12) and recalling that the video samples 
V- k are independent, yields 


Cm 


£ £ *(v t 

k 1 


V V 
ik' 


V 4 


a v, sv, 


d a t d 


(13) 


But for V- k exponentially distributed. 





( 14 ) 


20 



and the Information matrix becomes: 



* V * * 

" ¥ V nf 


a V, d V, 


d ^ 


(15) 


where N is the total number of pulses processed. Thus, the estimation 
accuracies can be determined from the mean return model given in 
Appendix A. Although analytical expressions can be developed for 
the partial derivatives and even for the information matrix for a 
particular set of parameters, the resulting equations give little 
insight to the behavior of the errors, nor do they significantly 
simplify the calculation of B. Therefore, both the partial derivatives 
and the computation of the covariance matrix are evaluated numerically. 
Examples for the GEOS altimeter parameter estimates are given in 
Appendix B. 


Reference 


[1] S. S. Wolks, "Mathematical Statistics," John Wiley & Sons, New 
York, 1963. 
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APPENDIX A 
MEAN POWER RETURN 


In this appendix, the radar equation for a distributed target 
Is derived and applied to the particular case of a satellite altimeter. 
The target Is modeled as a continuous distribution of Independent 
scattering points distributed In the three dimensions (t,v,u). Here, 
t Is range delay, v Is the Doppler shift, and u Is the sine of the 
angle to the scattering point measured from the antenna boreslght. 

Since the scattering points are assumed to be Incoherent, the total 
power Is simply the Integral over delay, Doppler shift, and angle 
of the differential power received from each point. The altimeter 
equation has been derived In various forms elsewhere (e.g., Barrlck 
[1,2], Harger [3] or Brown [4], This derivation Is Included for 
completeness, and to define the notation used to obtain an approximate, 
closed-form solution for the shape of the mean power return. 


Theoretical Developmen t 


From the theory of high resolution radar (e.g., Rlhaczek [5] 
nr Deley [6], the power received from a mismatched filter can be 
Mitten in terms of the cross-ambiguity function. That Is, let the 
received narrowband signal be represented In complex form as: 

S f (t) - /FT* g(t-T) exp[j2ll(f e -v)(t-T)] (A. 

where t Is the round-trip delay to the scatter, 

v Is the Doppler shift produced by the scatterer. 


f is the carrier frequency, 

g(t) is the complex modulation Impressed on the carrier, 
and E r Is the received signal energy. 


In the above expression it Is assumed that the complex modulation has 



been normalized to unit energy. 


If the receiver Is matched to a complex modulation w(c). a Doppler 
frequency $, and a delay, t, then the video power out of the filter 
at time t Is given by: 


V(t) • E r G q B | X gw (t-t. v)| 2 

where | X gw |‘ Is the normalized cross-ambiguity function defined by: 
x gw (itv) '/^ g(0 w*U+i)exp[-J2irvO dc 


(A. 2) 


(A. 3) 


and 

G 0 B • /|w(f)| 2 df - / |w(t}| 2 dt 


Is the product of receiver gain, times Its noise bandwidth. 

that has been normalized so that 
gw 




1 


Note 


If the scattering point has differential cross-section, do, and Is 
located at angular position u (In sine theta space) then the signal 
energy can be computed from the radar equation (Skolnlk [7]). That Is: 


m fjGA<u) do x 2 GA(u ) 
r 4^ 4AR 7 4* 


(A. 4) 


where 


Ej Is the transmitted energy, 

G Is the antenna gain on boreslght, 

A(u) is the one way loss factor when the target is not 

on boresight, 

X Is the R. F. wave length 


A-2 



do 1$ the differential cross-section at the scattering 
point, and Is generally a function of all three 
parameters (t, v, u), 

R ct/2 Is the range to the scatterer, 

and l Is the total system losses (except for processing 

mismatch, which Is Included In the cross-ambiguity 
function) . 


Combining Equations (A. 2) and (A. 4), gives the differential 
power received from a scattering point located at (t, v, u): 


dV(t) 


* T C X LC B .7 , 2 

— — r 2 *-1^2 | x ( T -t, V0)| d O (T»V,u) (A. 5) 

(An) 3 R* 


Integrating this equation over the three variables, and adding the 
receiver noise, yields the total average video power for a distributed 
target: 


v(t) 


7 7 

EC X LC B 
T o 

(4n> 3 


♦ kT oW 


f f I" I X gw ( T_t * V ' 0 ! 2 d0 (T » 

u *0 v : j ■-« R 


V,u 


(A. 6) 


where 

and 


k Is Bolzman's Constant, 

T q Is the reference temperature, 

F n Is the receiver noise figure. 


Equation (A. 6) is a fairly general representation of the power 
received from a distributed target. To apply It to the satellite altimeter, 
one must now compute an expression for the differential cross-section 
do(x, v, u) which Is appropriate to the altimeter. The geometry being 
considered Is shown in Figure A.l. 


Consider a scattering point located at height, h, above the 



But the second term Is a semldef Inlte quadratic form so that, 

dV<»bY 1 BD > o ( B . l9 ) 

and 

Q ( 0 n+l ) < Q(0 n > (B.20) 

Thus, the Iterative solution (eq. ( B . 1 6 ) ) converges to c local minimum 
for small enough k. Note that although • rather crude approximation 
was made to obtain eq. (B.16), the atnve discussion shows that the 
only criterion (for local convergence) Is that k be chosen small 
enough. Experience with the altimetry problem Indicates k on the 
order of .5 to .8 Is often times sufficiently small. 

APPROXIMATE VARIANCE OF THE ESTIMATE 


An approximate estimate of the covariance of the parameter 
estimates, 0, Is obtained by expanding eq. (B.5) about the true 
parameter value ». Then In analogy to equation (B.13), one has 

0 w 9 - [rc T (6)]" 1 e(0) (B.21 ) 


Where the above equation holds asymptotically for large numbers of 
samples. Again, assuming many samples, v r r (0) Is approximated by 
its mean value: 


Pe T (0) * EPe T (0) 

.T 


BB 


(B.22) 


where B is evaluated at the true parameter value, 0. Note, to obtain 
eq. (B.22), it has been assumed that 


Hence 


EVj - Vj (0) 

0 * 0 - (BB T ) _1 BD 


(B.23) 


(B.24) 


B-5 



(B.25) 


Again, assuming EVj“Vj , 0 Is (asymptotically) unbiased, and 

B 0 - 0 

The covariance matrix of the parameter estimates Is: 

* » 

Zq - E(0-0) (0-0) T 

T -1 T T -1 

- (BB ) B B*(BB ) 

For Independent samples, the covariance of D Is given by: 
5^ - Diag (o^ 2 /^ 2 ) 

where "Dlag" denotes a diagonal matrix. 

In the special case, of an MLE estimate for an exponential 
distribution, 

o Vj - Vj - ^ 

and equation (B.26) becomes: 

£ 0 - (bbV 1 

SUMMARY 

The important equations In the above discussion are: 
Iterative Solution 

0,-0 - kWD 

n+1 n 

W » (bb t ) -1 B 


B-6 


(B.26) 

(B.27) 

(B. 28) 
(B.29) 

(B. 16) 
(B. 30) 



B - (v'Vvj , Y l PV K ) (B.9) 

D- ((Vj-vp^ * l ... <V K -V K )^’ l ) T (B. 10) 


where for an exponential distribution maximum likelihood estimate: 



(B.7) 


but for a MMSE , 

Vj - Uj (B.8) 

where Uj is the weight applied to each sample residue squared. 

In the above expressions, B and D are evaluated at the current 
estimate, 0 n - 


Approximate Error Analysis 


E* - (bb T ) b2^ b T (bb T ) _1 

with 

^ - Diag (o v 2 /s> 2 ) 


Special Case 


For an exponential distribution, 


o = V 

V J J 


(B.26) 


(B.27) 


(B.28) 


and finally for an exponential distribution MLE , 
and 




T -1 
(BB T ) 




(B.29) 
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mean sea surface. Then, If o° (e) Is the mean cross-section per unit 
surface area (at angle e relative to radar) and If p h (h,o n ) Is the 
probability density function of the distribution of scatterers with 
height, (scaled to RMS wave height e (l ) then the differential cross- 
section at that point Is: 

do ■ o° (0) p^( h • o^) dV (A. 7) 

where dV Is the differential volume at the point of Interest. Now, 

In the spherical coordinate system (R.e,*) located at the satellite, 
the differential element of volume Is: 

dV « R sin e dR de d* (A. 8) 


The final step, then Is to compute t, v, u and h In terms of 
R, e and substitute the results Into Equation (A. 6) and then 
substitute Equation (A. 8) Into Equation (A. 6) to obtain the final 
expression. But from Figure A.l, it Is seen that 


2R 

c 

2V 


\ 

sin 0 


sin 0 cos iji 


and 

h ■ 

( (H Q -R cos e+a e 

, 2 ♦ R 2 sin 2 

where 

V 

is the satellite 

velocity, 


"0 

is the satellite 

altitude, 

and 

°e 

is the radius of 

the earth. 


(A. 9) 
(A. 10) 


Before making the indicated substitutions, the equations will be 
simplified by making some approx 'nations. That is, it is assumed 
that the beamwidth of the altimeter is narrow (less than 5°) so that 
small angle, and related consistent approximations can be made. 
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Figure A.l Geometry for Computing Mean Power Return 


Thus, It Is assumed that: 


u * sin e << 1 , 

R • H Q ♦ aR, 
aR « H Q 

aR << a„ 

e 

With these approximations, one can write 
fUft . m . T + 2A R 

r " e + “ V. c 
He 2 H 0 

h « -AR + -y" (1+ j 2 ) 


(A. 11) 


(A. 12) 


do « o°(0) p h (h.c^) H o 2 0 dAR dfl d« 

2Ve 

ano v — y— cos ^ 

Thus, the average video power can be written finally as: 


v(t) 


E G 2 A 2 Lfl°G B 
T o 

(4ll) 3 H 3 (1+H / a ) 
o o e 



dAR 



[H (1+H /a )6 dG] 

o o e 


Y(0)A ? (fl)p h (- AR 


H e 2 v 

+ "V <1 + H o / V» a h) 


/ 



, 2V9 
' X 


COS ijf - 



2 


+ kTFG B (A. 13) 

o n o ' 

where Y ( 0 ) is defined by o°(e) * o°y(e) and defines the variation of 
o°(e) with angle. 
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Ap proximate Evaluation 


A closed form solution for the mean power return can now be 
obtained from the triple Integral of E^-ation (A. 13). To achieve 
such a solution requires some additional approximations and simplifying 
assumptions. First assume that the cross ambiguity function Is 
essentially constant over the doppler variation of the scattering 
points, then 



Further assume that the variation of the cross ambiguity 
function with range delay can be approximated by a Gaussian function, 
thus 



n(t L , o T ) 


(A. 15) 


where n( • »o ^ ) Is the normal density function and * RMS pulse width. 


Now from Barrick's [1] mooel , the wave height density is also 
normal , hence 


p h (h, V “ 

Letting 

H 

H - H (1 + T 2 ) 

o a 



(A. 16) 


(A. 17) 
(A. 18) 
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dx ■= HOdP 



Substituting Equations (A. 14) through (A. 18) Into the triple Integral 
of Equation (A. 13), denoted Ij, yields. 


I 


3 




dxY(x)* 2 < x > 



dAR 


(A. 19) 


H(x - AR, 



<t-T 0 )c/2, o R ) 


where 


°R 


CO 

■ RMS range resolution. 


Thus the Integral over aR is simply the convolution of two normal 
densities having different means and variances. As may be shown, 
this convolution produces a normal density with a mean equal to 
the sum of the means, and a variance equal to the sum of the 
variances. Letting 

(A. 20) 


(A. 21) 
(A. 22) 

Now assuming that the two-way antenna pattern function can be 
approximated by a Gaussian function 


y ■ 3}Ct « range at time t, 

and y “ 1&CT = range to the mean sea surface 

o o 

and performing the indicated convolution of Equation A. 19 gives 

HTT 2 

3/2 /»® o 

X 3 ■ (2tt) Or j y < x ) a (x) TjCy - y Q - x, o £ ) 

0 

where o, 2 - + 


2 

A (0) ~ exp[-(4 In 2) (0/ 6 BW2 ) 2 3 
where ®BW2 ~ *- wo ' wa y 3 dB beamwidth. 


(A. 23) 
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and from Barrlck's [1] model, the surface shaping function Is 


Y<0) ~exp[-(4 ln2)<e/e BW2 ) Z ] 


(A. 24) 


wnere 0^ " the 3 dB spread of the slope distribution. 

Then, utilizing the change of variable given by Equation A. 18, the 
product of the surface shaping function and beam pattern can be 
written as 


Y(x)A (x) - exp(-nx) 


(A. 25) 


where 




3^ 

dW2 


e: - e 


- 1/2 


(A. 26) 


Is the effective 3 dB beamwidth, and 


ji - (8 1n2)/(H0 e ) 


(A. 27) 


is the exponential decay factor parameter. 


The remaining Integral of Equation A. 21 can now be expressed as 

.2 


HTT - 


r 8 2 

F(z) - J Y(x)A (x) T|(z - x, o e ) dx 


(A. 28) 


V2n 


h~ / cxp b(^r) I d * 


where it has been assumed that HTT /8» o - The remaining integral can 

e 

be evaluated by completing the square in the exponent. This yields: 
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Tiz) 


•x P r%(0 # H> 2 -*$»l^ 



0 


dx (A. 29) 


So that finally 


F(x) - cxp[%(o # n) 2 -*n] * <o e ’ V> {A 30) 

where <t> Is the cumulative normal distribution function. Thus, 

F(z) describes the shape of the mean power return as a function of 
range measured from the mean sea surface. 


The final expression for the mean return may be obtained by 
substituting (A. 30) Into (A. 21) and the result Into (A. 13). Before 
doing this, note that the peak slgnal-to-nolse ratio at the output 
of the receiver Is given by: 


(2tt) 3/2 e t G 2 X 2 Lo^o° 

( 1 + £“) k V„ 

e 

Thus, the mean return may be written: 

V(y) - k T FGB (a F (y - y ) + l) 
o n o \ o / 


where the function F(.) Is given by Equation (A. 30). Note that If 
the noise power Is known (XT 0 F n G 0 B), then the mean return depends 
explicitly on the four parameters: 


a - signal-to-noise ratio 
y Q - altitude to the mean sea surface 
o g - "effective" wave height 
and w - exponential decay factor. 

The derivation of Equation (A. 32) did not Include the effects of 

antenna pointing error. These effects may be Included by noting 

2 

that pointing error causes the two-way beam pattern, A to depend 
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both on the off-nadir angle, 0, and the azimuth angle In particular, 
for an off-nadir angle < we have: 

9 # 2 

A (0,t) - exp [-(A In 2)(0/e BW2 ) ] (A. 33) 

where 0 Is the angle between the antenna boreslght and a line to 
the surface scattering point located at angular position [8,+). 

0 Is determined by the equation: 

Cos p - Cos 0 Cos? + Sin0 Sin§ Cosf (A. 34) 

Making small angle approximations yields 

• 0 2 „ e 2 + S 2 ♦ Cos* (A . 35 ) 

Since A2 now depends on the azimuth angle It must be factored 
Into the evaluation of the Integral over g* In Equation A. 14. 

Thus, It Is necessary to evaluate: 

•4 (A. 36) 

' t.^l 2 A 2 (6) 


where 



and I Q Is the modified Bessel function of order zero. 

Equation (A. 37) may be further simplified by noting that for 
arguments less than 1, I may be approximated to about 3 percent accuracy 
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by: 


I o (t) - e 



(A. 38) 


Then, 



(0) « exp 


(4 1n2)g 2 
0 BW2 



where 



(A. 39) 


The development proceeds Identically now, except that A c Is used 
everywhere A^ appeared previously. Thus, the mean return model 
given In Equation (A. 32) Includes pointing error If the slgnal-to- 
nolse ratio Is modified as: 


a' - * cxp[ - (A In 2)(^— ) 2 ] /« . 

. 0 BW2 * 

i 

and the modified two-way beamwidth In Equation (A. 39) Is used In 
the definition of the exponential decay factor (A. 27). 
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APPENDIX B 

MATRIX FORMULATION OF THE MLE AND MMSE ALGORITHMS 
INTRODUCTION 


In this Appendix, the Iterative equations for obtaining the 
maximum likelihood estimator (MLE) for an exponential distribution 
and the minimum mean square error (MMSE) estimator expressions for 
the covariance matrix of the estimators Is given. 

DEVELOPMENT 


As was done In earlier studies (e.g., ref. [1], It Is assumed 
that K Independent data samples of a distributed target are taken 
at the output of a square law detector. Denote these samples as: 

V^, j - 1 ••• K 

Further, assume that a parametric form of the mean return model Is 
known. That Is, 


v j(6), J - 1 K 

where 0 Is a vector of parameter values to be estimated. 

Under these assumptions, the logarithm of the likelihood function 
Is given by: 

A(0) - - E tin Vj + (V^)] (B.1 ) 

1-1 

The MLE estimate of the parameters, are those values, 6 , which 
maximize eq. (B.l). Equivalently, let the MLE penalty function be: 

Q (0) “ -A(0) 

MLE 

K 

- £ [in V + (V./V )] (B.2) 

J-l J J J 

A 

and 0 minimizes Q^f* 

I 
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In a similar manner. If the quadratic penalty function Is 
defined as: 




* _ 2 

E <<V v .)/u.r 

J-l J 3 1 


(B. 3) 


Then the (weighted) minimum mean square error parameter estimates 

. n 

are those which minimize (where the weights are Uj , J-1....K) 


If Vis used to symbolize the gradient (w.r.t.e ) operator, then, 


V 



(B.4) 


A necessary condition on the estimators Is that the following M 
.'tlons be satisfied simultaneously. 

*(0) - 7Q<0) 


(B.5) 



Note that by differentiating eq. (B.2) and eq. (B.3), one may write a 
unified version of eq. (B.5) as: 


c(6) - E [(V -V.)/(v 4 Z )] ^ v 


J-l 


J J J 


’J 


(B.6) 


where to obtain a MLE estimate, one lets 



and to obtain a MMSE estimate, one lets 

V j ’ “j (B.8) 

The matrix form of eo. (B.5) Is obtained by letting B be the K 
by M matrix defined as 

> - tv,' 1 ?*,. •••V' v,) 


(B.9) 


and letting D be the K dimensional norma'lzed deviation vector, 


(B.10) 


Then the necessary condition becomes 

e (6) - BD - 0 (B.ll ) 

To find an iterative solution to eq. (B.5) or (B.ll), let # n be 
the current approxiniation, and expand eq. (B.5) in a Taylor series about 
« n , thus, 

0 « c ( 0) « e(0 n ) + [7t T (0 n )3 (O-0 n ) (B.12) 
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and 


e " o - [Pc T <e )]’ 1 «(e ) 

n n n 


(B.13) 


Differentiating eq. (B.6), one finds: 


Assuming reasonably good starting values, one has V j - V ^ 
second term Is small compared to the first so that; 


- 0 


yc 




T 

(BB T ) 


Thus the Iterative equation becomes 


Vi ■ V k - <“V »0 


In the above equation, B and D are evaluated at the current 
estimate 0 . The constant k Is chosen less than unity. In 
for k small enough the above procedure converges to a local 
That Is, 


«Vi> ' Q( V + [rT Q(6 » )3 <6 n+rV 


or 


+ higher order terms 

Q<e M l > ” «<V + ,T<e »> 

- Q(0 n ) - kD T B T (BB T )" BD 


(B. 14 • 

, and the 


(B.15) 


(B. 16) 


fact, 

minimum. 


(B. 1 7) 


(B.18) 
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